The Extension of Rod-Coil Multiblock Copolymers 
and the Effect of the Helix-Coil Transition 
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The extension elasticity of rod-coil mutliblock copolymers is analyzed for two experimentally acces- 
sible situations. In the quenched case, when the architecture is fixed by the synthesis, the force law 
is distinguished by a sharp change in the slope. In the annealed case, where interconversion between 
rod and coil states is possible, the resulting force law is sigmoid with a pronounced plateau. This last 
case is realized, for example, when homopolypeptides capable of undergoing a helix-coil transition 
are extended from a coil state. Both scenarios are relevant to the analysis and design of experiments 
involving single molecule mechanical measurements of biopolymers and synthetic macromolecules. 



PACS numbers: 61.25.Hq, 61.41. +e, 87.15. He 

With the advent of single molecule mechanical mea- 
surements it became possible to study the force laws char- 
acterizing the extension of individual macromolecules [Q . 
In turn, these provide a probe of internal degrees of 
freedom associated with intrachain self assembly or 
with monomers that can assume different conformational 
states. A molecular interpretation of the force laws thus 
obtained requires appropriate theoretical models allow- 
ing for the distinctive "internal states" of each system. 
The formulation of such models is a challenging task in 
view of the complexity and diversity of the systems in- 
vestigated. These include DNA the muscle protein 
titin |3) and the extracellular matrix protein tenascin [Q , 
the polyscharides dextran Q and xanthan [01 as well as 
the synthetic polymer poly(ethelene-glycol) 0. 

In this letter we consider two unexplored yet accessible 
systems where a detailed confrontation between theory 
and experiment is possible. In particular, we present a 
theory for the extension force law of multiblock copoly- 
mers consisting of alternating rod and coil blocks. Two 
scenarios are considered, both focusing on the equilib- 
rium force law of long chains undergoing quasistatic ex- 
tension. In one, the architecture is "quenched" that is, 
the block structure is set by the chemistry and no inter- 
conversion is possible. Such is the case, for example, for 
segmented polyurethans ||. In the second scenario the 
monomers can interconvert between coil and rod states. 
The "annealed" architecture is realized, for instance, in 
homopolypeptides capable of undergoing a cooperative 
helix-coil transition p)-|Tl|]. In this system, the highly 
rigid helical domains play the part of the rod blocks. As 
we shall see, the two scenarios lead to distinctive force 
laws (Figure). The force law of the quenched case is 
characterized by an abrupt change of slope. This arises 
because the rod blocks are more susceptible to orienta- 
tion by the applied tension. In the annealed scenario, the 
extension of a chain that is initially in a coil state leads to 
a sigmoid force profile exhibiting a pronounced plateau. 
The plateau is traceable to a one dimensional coexistence 
of helical and coil domains where the chain extension fa- 
vors the helical state because of its low configurational 



entropy. These results suggest that force measurements 
can be used to probe the block structure of polymers with 
quenched architecture. In the annealed case the force law 
provides a direct measure of the thermodynamic param- 
eters involved. From the view point of the growing field 
of single molecule mechanical measurements, these re- 
sults are helpful in exploring the diagnostic potential of 
these techniques. The discussion is of interest from the 
polymer physics perspective because the theory of poly- 
mer elasticity focuses on the case of flexible homopoly- 
mers modeled as random walks, or self avoiding random 
walks, with a constant step length and no internal de- 
grees of freedom jllj]. In this context, rod-coil multiblock 
copolymers with quenched architecture may be viewed as 
a special case of heteropolymers incorporating monomers 
of different sizes. The analysis of the annealed case sup- 
plements the isolated discussions of the effect of internal 
degrees of freedom on the extension elasticity |lj [ll| . 

We first consider the elastic free energy F e i and the 
tension / = 8F e i/dR of a quenched rod-coil multiblock 
copolymer with an imposed end-to-end distance, R. The 
chain incorporates N monomers of identical length a that 
form yN rod blocks and yN coil blocks such that the 
number of monomers in the rods is ON . For simplicity we 
assume that the coil and rod blocks are monodispersed, 
and that the number of monomers in a rod block is thus 
0/y. We focus on the case of long rod blocks, 0/y ^> 1. 
The distinctive feature of rod-coil multiblock copolymers 
is that the rodlike blocks, the longer "monomers" , are 
oriented by the applied tension before the alignment of 
the shorter monomers becomes significant. The freely 
jointed chain model, the macromolecular analog of the 
Langevin's theory of paramagnetism, enables a quanti- 
tative description of this situation jl^Jlfj]. The length 
of the monomer is the counterpart of the magnetic mo- 
ment, the applied tension is the analog of the magnetic 
field while the end-to-end distance of the chain corre- 
sponds to the magnetization. A monomer of length I is 
assigned an orientational energy of —flcoscj) where ip is 
the angle made by 1 with respect to f. For a flexible 
homopolymer the end-to-end distance is R = NaL(x) 
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where L(x) = cotha; — x^ 1 is the Langevin function and 
x = fa/kT. In our situation the multiblock copolymer 
is viewed as a collection of two types of non-interacting 
"dipolcs": (1 — 9)N monomers of length a and yN 
"monomers" of length 6a/ y. The reduced end-to-end dis- 
tance, r = R/Na is thus 



fa/kT 




FIG. 1. The force laws characterizing the extension of a 
rod-coil multiblock copolymer with a quenched (a) and an 
annealed (b) architectures. The parameters for (a) are 6 — 0.2 
and y = 10~ 4 . The continuous line depicts the exact force law 
while the dashed line describes the approximate expression 
discussed in the text. The chain configurations in regions (I) 
and (II) are schematically depicted above. The parameters 
for (b) are s = 0.7 (dashed line) and s = 0.8 (continuous 
line). In both cases a = 10~ 4 . The dotted line depicts the 
S m i X — approximation for s = 0.8. f co , r^, ri/ 2 and rjj 
correspond to s — 0.8 case. 

r = (l-6)L(x) + 6L(6x/y). (1) 

Three regimes are involved. When i < 1 and Ox/y <C 1 
both the rods and the a monomers are weakly ori- 
ented. Since L(z) z/3 for z < 1, we find fa/kT ss 
3r/(l - 9 + 9 2 /y) and F e i/NkT w 3r 2 /2(l - + 6 2 /y). 
For stronger extensions x <C 1 while Ox/y 3> 1 in 
other words, the a monomers are only weakly aligned 
but the rods are almost fully oriented with /. Focusing 



on the case of long rods, y <C 1, when L{9x/y) » 1, 
we obtain fa/kT » 3(r - 0)/(l - 0) and F el /NkT « 
3(r - 0) 2 /2(l - 0) + 3y/2. The first term in F e , reflects 
the Gaussian stretching penalty of the (1 — 9)N "short" 
monomers. It allows for the full alignment of the rods 
and the resulting modification of the imposed end-to-end 
distance experienced by the flexible blocks. The second 
term reflects the Gaussian penalty associated with the 
full alignment of the yN rod blocks. The crossover be- 
tween the two regimes occurs at r co 9 + y(l — 9)/ 9. 
Eventually, upon further extension, the a monomers also 
approach saturation, i» 1 and 9x/y ^> 1. In this last 
regime fa/kT w (1 - + y)/(l - r) and F el /NkT » 
(1-6 + y) ln[2(l - 0)/3(l - r)\ + 3y/2 + (1 - 0)/6. Since 
the approximate expressions for / do not match at the 
boundary between the two last regimes, the crossover 
between them, f^,, may be determined by the require- 
ment that fa/kT « 3(r - 0)/(l - 0) » 1 leading to 
»Q, s» (1 + 20) /3. Altogether, the plot of fa/kT vs. r 
for a rod-coil multiblock copolymer is distinguished by 
an abrupt change of slope at r co . In turn, this feature 
provides a useful diagnostic for the architecture of the 
polymer. 

In the annealed scenario, monomers may interconvert 
between the rod and coil states. Accordingly, 9 and y are 
no longer constants set by the chemistry. Rather, their 
values vary with the temperature T and with R. 9 and y 
are determined by minimizing the free energy per chain 
Fchain with respect to y and 9 for a given R. F c h a in al- 
lows, in addition to F e i , for the mixing free energy F m i X 
of the one dimensional mixture, along the chain trajec- 
tory, of the rod and coil blocks. Our discussion focuses 
on polymers undergoing a helix-coil transition. It is di- 
rected at the coupling of the helix-coil transition with 
the extension for the simplest possible situation, that is 
a homopolypeptide capable of forming an a- helix § 0. 
This form of intrachain self assembly is due to intrachain 
hydrogen bonds. The i-th monomer, residue, forms H- 
bonds with the (i + 3) and (i — 3) monomers. Altogether, 
a helical domain consisting of n monomers involves n — 2 
H-bonds. The persistence length of an a-helix is of order 
of 200nm and a helical domain may thus be considered 
as a rod block. It is convenient to discuss this system 
in terms of "bonds" between two adjacent monomers. 
For simplicity we retain the approximation setting the 
monomer size in the rod and coil state equal to a. Us- 
ing the coil state as a reference, each helical bond is as- 
signed a free energy A/. It reflects the contribution of 
the intrachain H-bonds, the change in solvation due to 
the formation of intrachain H-bonds and the associated 
loss of configurational entropy. The helix-coil transition 
occurs at the transition temperature T* when A/ = 0. 
Above T*, A/ > while below T*, Af < 0. Terminal 
helical bonds, at the boundary of a helical domain, incur 
an additional free energy penalty Af t > since the ter- 
minal monomers lose their configurational entropy but 
do not contribute H-bonds. A/ t plays the role of inter- 
facial free energy associated with the helix-coil domain 
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boundary. Traditionally, the theory of the helix-coil tran- 
sition is formulated in terms of the Bragg-Zimm param- 
eters s = exp(— Af /kT) and a = exp(— 2A/ t /fcT) where 
a « 1CP 3 — 10~ 4 varies with the identity of the residues 
but is independent of T. The helix-coil transition gives 
rise to a sigmoid 9 vs. T plot with a characteristic width 
~ X^er 1 / 2 . The usual discussion of the helix-coil transi- 
tion is based on transfer matrix methods. For our pur- 
poses it is convenient to recast it in terms of the free 
energy of the unperturbed chain, F® hain . For long chain, 
when the N — > oo limit is applicable, F® hain reflects 
three contributions. Each of the 9N helical bonds con- 
tributes Af while every one of the 2yN terminal bonds 
contributes A/ t . These two terms are supplemented by 
the mixing entropy, S m i X , associated with the different 
possible combinations of 9N helical bonds and 2yN ter- 
minal bonds. The number of possible configurations, W 
is WtbWtc where W t h = {J^) is the number of ways of 
placing yN terminal bonds among 9N helical bonds and 
Wtc — i^ 1 ~ y N N ) is the statistical weight associated with 
the placement of yN terminal bonds among (1 — 6)N 
coil bonds. S m i X = kh\W can be expressed as S m i X = 
eS mix (y/9) + (l-9)S mix (y/(l-9)) where ±S mix (X) = 
-XlnX - (1 - X)ln(l - X). Here j^S mix (y/9) is the 
mixing entropy of the terminal monomers among the rod 
monomers while j^S m i X (y / '(1 — 9)) is the mixing entropy 
of the terminal monomers among the coil monomers. Al- 
together 
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The equilibrium conditions 8F® hain /dy = and 
dF a hain /d9 = lead respectively to y 2 = (9 - y)(l - 
9 - y)a and to (1 - 9)(9 - y) = 9(1 - 9- y)s. In turn, 
these yield the familiar results for the unperturbed chain 

and, in particular, B=\ + \ \J ~jJ^^ s- 

The full analysis of the coupling of the helix-coil tran- 
sition with the extension of the chain involves an aug- 
mented free energy per chain, allowing for the extension 
penalty F cha in = F° ham + F e i @. An analytic solution 
yielding the qualitative features of the force law, as well 
as the important length and force scales, is possible by 
setting S m i X = Within the S m i X — approxi- 

mation the interfacial penalty imposes coalescence of the 
helical domains leading to the formation of a helix-coil 
diblock copolymer (y = 0). As a result, the helix-coil 
transition takes place as a first order phase transition. 
The corresponding free energy is 
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The equilibrium condition dF c h a in/ 89 = for a given r 
yields 9 = 1 - (1 - r)/y/l - 2A//3fcT. In turn, this de- 
fines a critical extension, = 1 — a/1 — 2A//3fcT such 



that for r < 7"l, — and for r > r^, 9 is specified by 
9 = (r — tl)/{\ — r£). The equilibrium free energy of the 
chain for r < tl is F c h a in/NkT w 3r 2 /2 while for r > tl 
it is F chain /NkT « 3r L (r - 1) + Af/kT. This regime 
lasts until 9 = 1 is attained at rjj — 1. Within this sim- 
ple minded picture, further extension is impossible and 
for rjj > 1, F c hain = oo. The corresponding force law in- 
volves three regimes: (i) A linear response regime, where 
the Gaussian elasticity is operative and fa/ NkT w 3r, 
occurs while r < rz, and 9 = 0. (ii) A plateau with 
f co a/NkT « 3r\L occurs in the range < r < rjj when 
the coexistence between the helix and coil states follows 
the lever rule 9/(1 — 9) = (r — ri/)/(rxj — r). The mid- 
point of the plateau, corresponding to 9 — 1/2, occurs 
at r x / 2 = (ru + tl)/2. (hi) A steep increase in force oc- 
curs at r = rjj, when the assumed nonextensibility of the 
fully helical chain comes into play. Note that in reality 
the helical domains are not perfect rods. High applied 
tension may cause extension by breaking the H-bonds. 
Within the S m i X — approximation the plateau corre- 
sponds to a first order phase transition involving the co- 
existence of a helical phase and a coil phase. Such phase 
transition is prohibited in one dimensional system expe- 
riencing short ranged interactions ]l8j . When one allows 
for S m i x > the plateau in the force law exhibits a weak 
increase with r, instead of the r° dependence predicted 
by the S m i X — approximation. Yet, as we shall see, 
the center of the plateau occurs at 7*1/2 and its height at 
this point is f co a/kT « 3r£. Furthermore, the qualita- 
tive Af dependence of both the height and the width of 
the plateau is retained. In particular, the width of the 
plateau increases while its height decreases as Af — > 0, 
or equivalently, as s — ► 1. Since at the vicinity of the 
transition temperature s ~ (T — T*)/T*, these scenarios 
may be explored by change of T. 

The complete analysis of the problem in terms 
F c ha m = F° hain + F el involves different regimes dis- 
tinguished by the applicable form of F e i, as discussed 
in the first part of this letter. For weak extensions 
F el /NkT « 3r 2 /2(l - 9 + 9 2 /y) and the equilibrium 
conditions dF c h a in/d9 = and dF c h a in/dy — lead to 
dF° cham /d9 = -8F el /d9 and dF° ham /dy '= -8F el /dy. 
In turn, this leads to y 2 = (9 — y)(l — 9 — y)oK a 
and to (1 - 9)(9 - y) = 9(1 - 9 - y)sK s where K a = 
exp[(dF el /dy)/NkT] and K s = exp[(dF el /89)/NkT]. 
(8F e i/d9)/NkT is the increment in the elastic free en- 
ergy per monomer associated with adding one residue 
to a helical block for a constant y and (dF e i/dy)/NkT 
is the price of creating an extra helical sequence while 
maintaining a constant 9. For weak extensions, when the 
elastic penalty is a weak perturbation, K a and K s may be 
approximated by K a = exp[-r 2 (3 2 /y 2 (l - 9 + 2 /y ) 2 } 
and K s = exp[r 2 (2V2yo - 1)/(1 - #0 + 2 /y Q ) 2 } where 
9q and yo correspond to the unperturbed chain. K a 
is a decreasing function of r while K s is an increas- 
ing function of r. Thus, the extension favors the tran- 
sition by increasing s and enhances the cooperatively 
by decreasing a. This analysis demonstrates that the 
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Smix = approximation overestimates F e i and thus /. 
fa/NkT » 3r/(l — 9 + 9 2 /y ) provides abetter estimate 
for /. However, in accordance with the Le Chatellier 
principle |p| , the actual force law is even weaker since 
9 increases with r while y/9 decreases. For stronger de- 
formations, when r > r co the appropriate elastic free 
energy is F e i /NkT « 3(r - 6>) 2 /2(l - 6) + 3y/2. In this 
range it is possible to obtain the force law at the vicinity 
of 9 = 1/2. The equilibrium condition dF c h a in/dy = 
leads to y 2 = (9-y)(l-9-y)a/e 3 / 2 . Since a <C 1 this re- 
duces, in the vicinity of 6 = 1/2, to y ex \J 9(1 — 9)a/e 3 / 2 . 
Utilizing this relationship, and following expansion of 
the logarithmic terms we obtain F c h a in/kT ps 3(r — 
(9)72(1 - 6) + 9Af/kT - 2^/9(1 - 9)a/e 3 / 2 . In turn, 
dF c hain/d9 = recovers the expression for 9 as obtained 
within the S m i X = approximation, 9 » (r— rj,)/(l— Tl)- 
This allows us to express -F^/iam & s a function of r and to 
obtain / in the vicinity of 9 = 1/2 

fa f co a a 1 ^ 2 (r ^ r 1/2 )(l ~ tl)- 1 

kT ~ kT e 3/4 ^{r-r L ){ ru -r) ' [ ' 

The force at r±/2 is f co as suggested by the S m ix — 
approximation. However, when the full F c h a in is allowed 
for the / ~ r° behavior of a perfect plateau is lost. The 
second term in (Q) ensures that / increases with r. On 
the other hand, since a « 10 -4 , the a 1 / 2 prefactor in this 
term imposes slow variation of / in the vicinity of ri/ 2 - 

Our discussion of the elasticity of rod-coil multiblock 
copolymers reveals two scenarios. In the quenched case, 
the force law exhibits a sharp change of slope while in the 
annealed case the force law is sigmoid. It is important 
to note that our analysis focused on the thermodynamic 
force law as obtained for long chains and quasistatic ex- 
tension. Effects due to the finite length of the chains 
and the finite rate of deformation may modify the ex- 
perimentally observed force curves jl5). The analysis of 
the coupling between the helix-coil transition and the 
chain extension focused on the case of homopolypeptides 
capable of forming a single strand a-helix. The results 
can be however applied, with some modifications, to a 
much wider class of systems. It is relevant, for exam- 
ple, to the interpretation of the sigmoid force curve ob- 
served for poly(ethylene-glycol) which was attributed to 
the formation of water- mediated single stranded helix . 
Furthermore, the analysis can be extended to the case 
of polymers forming multistranded helices such as colla- 
gen. In this context it is of interest to note that sigmoid 
force curves were indeed observed in early experiments 
involving fibers formed by fibrous proteins, collagen, ker- 
atin etc., that form helical structures [^9|. Finally, one 
should note the recent development of facile synthesis 
route for homopolypeptides and block copolypeptides of 
well defined architecture pof . This provides a convenient 
method for synthesizing rod-coil multiblock copolymers 
with either quenched or annealed architecture. 
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